Mesoscopic magnetoelectric effect in chaotic quantum dots 
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The magnitude of the inverse Faraday effect (IFE), a static magnetization due to an ac electric 
field, can be strongly increased in a mesoscopic sample, sensitive to time-reversal symmetry (TRS) 
breaking. Random rectification of ac voltages leads to a magnetization fiux, which can be detected by 
an asymmetry of Hall resistances in a multi-terminal setup. In the absence of applied magnetic field 
through a chaotic quantum dot the IFE scale, quadratic in voltage, is found as an analytic function 
of the ac frequency, screening, and coupling to the contacts and floating probes, and numerically 
it does not show any effect of spin-orbit interaction. Our results qualitatively agree with a recent 
experiment on TRS-breaking in a six-terminal Hall cross. 
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In pursuit of effects that combine spin and charge, 
great attention is paid to mutual effects of magnetic and 
electric fields. Their magnetoelectric manipulation is in- 
teresting not only scientifically, but also for possible ap- 
plications. ^ Obviously, any magnetooptic or magneto- 
electric effect has its counterpart. For example, a fer- 
romagnet polarizes the spins of electrons and conversely, 
a current can exert a torque on a magnetization vector 
and switch a magnetic domain in random access mem- 
ory. ^ Similarly, the Faraday effect, a rotation of the po- 
larization of light by a magnetic field, has its inverse: a 
medium is magnetized by a beam of circularly polarized 
radiation. ^ We consider the unusual properties of this 
magnetoelectric effect in small non-magnetic samples. 

Classically, if a medium has a spatial inversion or time- 
reversal symmetry (TRS), the expansion of its free energy 
F does not have terms linear in E or H, respectively. 
Magnetoelectric effects appear only due to mixing of H 
and E, either as an HiEj term if both symmetries are 
broken in equilibrium (magnetized anisotropic medium) 

or in a higher order " HEE term" ^ HiEjE^^*'' if a strong 
electric field breaks TRS. In particular, this HEE term 
leads to magnetization AI — dF/dH, detectable in non- 
linear effects like second-harmonic generation or rectifica- 
tion^. The latter, known as inverse Faraday effect (IFE), 
is a static magnetization induced by perturbations at the 
frequency uj. "^'^ An electrical ac voltage, Kj, generated 
across a non-absorbing diffusive medium with mean-free 
time T induces an asymmetry of the dielectric tensor. 
Since it is also linear in magnetic field, the classical IFE 
magnetization flux, (pci, can be estimated,"^ 

(fid ^ a[s<S>o{eVu,)^mvl/{hujfc, ujt > 1, (1) 

where afg is the fine structure constant, ccfs ~ 1/137, and 
$0 = h/e is flux quantum. Importantly, this estimate 
could as well be obtained from the Joule heating and 
the asymmetry of the conductivity tensor. Since Eq. (1) 
contains small vp/c and the fine-structure constant, one 
naturally asks: can we enhance the magnetic response to 
the i?-field, which eventually breaks TRS? 

To show that it is possible to exceed Eq. (1), we propose 
to use the sensitivity of electronic transport through a 



mesoscopic sample to the broken TRS (Ref 7) as a detec- 
tor of magnetic flux. Indeed, disorder inevitably breaks 
a spatial symmetry in such a sample, and mesoscopic 
(sample-to-sample) fluctuations of transport occur on a 
fiux scale $ ^ <I>o- TRS breaking could be induced either 
by a flux $ of applied magnetic held, or by an IFE flux 
(fi created by additional ac perturbations. At $ = 0, an 
indirect transport detection of ip becomes possible in a 
multi-terminal setup with a separate pair of current- and 
voltage-probes^ measuring the Hall (or non-local) resis- 
tance Rh- Previously, Edelstein considered Rh oc as 
a signature of IFE in a noncentrosymmetric diffusive two- 
dimensional material with the mirror symmetry broken 
by spin-orbit interaction (SOI). ^ Due to large spin-orbit 
scattering time, Tso ^ t ^ l/o;, this effect is small, 
fso ~ ^o{eVui / Sf)"^ (h/ ujeft^q) ■ In contrast, here we con- 
sider a large fluctuational effect. 

Importantly, Rh 7^ even in small chaotic samples at 
zero field, but its random response to can be used as 
a sample-specific gauge for the fiux. Biittiker showed^" 
that the Onsager symmetry relations hold in mesoscop- 
ics, Rh{^) = RJj{—^)) where stands for the measure- 
ment with current- and voltage-probes exchanged; this 
prediction was confirmed in many linear transport exper- 
iments. 11^13 Only recently did Chepelianskii and Bouch- 
iat show that these relations are violated, when TRS is 
broken by additional ac perturbations. Since both ap- 
plied $ and induced lead to the asymmetry Rh 7^ i?^, 
data in Rh{0) — Rji{0) suggest a shift of zero magnetic 
flux, i.e. an effective IFE flux f 0. 

In this Rapid Communication we develop a theory that 
evaluates IFE in a multi-terminal chaotic quantum dot 
subjected to external ac perturbations at the frequency uj 
in the absence of applied magnetic field. First, we intro- 
duce a (sample-specific) gauge for magnetic fiux, using 
the response of Rh to small Then we find the TRS- 
breaking in Rh for the perturbed dot and evaluate the 
scale of fiuctuations of the induced IFE flux tf through 
its area: 

(fi - ±$c(eK;/£)^ eK; < e = max {huj, h/vd, T}, (2) 
where $c ^ 'i'o is the flux that completely breaks TRS, 
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Td is a typical dwell time of electrons in the dot, and T 
is the temperature. Equation (2), the main result of our 
work, does not have definite sign, which is typical for 
quantum effects. Similarly to the equilibrium persistent 
current or magnetic response of quantum dots, ^^'^^ it 
does not contain afs or large ep. Experimentally, a sam- 
ple can be magnetized not only by the flux <I> in equi- 
librium, but also by ac voltages at $ = 0. Specific de- 
tails of the sample rectify these voltages and lead to a 
random TRS breaking interpreted as an IFE magneti- 
zation. We expect this quantum interference effect to 
be measurable in any coherent sample of reduced dimen- 
sions. First we explain the model and major steps in the 
derivation of Eq. (2). Then we discuss how IFE is af- 
fected by screening, spin-orbit, and floating probes, and 
qualitatively compare our predictions to experiment. 

Model. We consider electronic transport through a 
multi-terminal A'^-channel chaotic quantum dot at the 
temperature T, see Fig. 1(a), and for simplicity assume 
the dot to be circular with radius L. Chaos is set due 
to either the diffusive motion of electrons, I = tvf <C i, 
or their random scatterings from the boundaries, L <^ I, 
so that the ergodic time ~ max {1,L/1}L/vf needed to 
explore the dot is short. ^^'^^ Therefore, electronic fo- 
cusing and direct trajectories, present in some ballistic 
structures, are absent. Each of four contacts has 
Na ^ 1 ballistic orbital channels and is characterized 
by Ua = Na/N. Additional coupling (not necessarily 
ballistic) of floating probes ^^^^"'^ is discussed in the end. 
ac voltages Va cos{ujt + 4>a) applied at the same frequency 
uj, but generally out of phase, are for convenience speci- 
fied by the vectors Va = V^e*"^" in a complex plane [see 
Fig. 1(b)]. A top gate with capacitance C is biased by 
Vo. 

Screening in the dot is accounted for by the uniform 
time-dependent potential U, and its higher spatial har- 
monics are suppressed due to short ergodic time. This 
potential is found self-consistently from charge conser- 




FIG. 1: (a) A gated chaotic dot in a four-terminal set-up, 
where the current through source-drain contacts 1 and 2 and 
the voltage drop between 3 and 4 yield Rh = -Ri2,34. IFE due 
to additional ac voltages Va cos{ijjt + (j>a) induces a random 
magnetization flux (p, detectable in Rh — RJj; (b) Biases are 
specified by Va ~ Vae"*° in the complex plane. IFE depends 
on the potential of the dot, [/ « Vb (or U') for weak (or 
strong) interactions. 



vation and gauge invariance. The contributions lead- 
ing to Coulomb blockade are small in 1/A^ ^ 1 and 
not taken into account. ^■^ If A = 2h^/mL^ is the mean 
level spacing, the ratio e^/CA defines the limits of 
weak/strong interaction. Using the dwell time of elec- 
trons, Td = h/NA, we normalize the flux of applied mag- 
netic field by the crossover flux that completely breaks 
TRS, = where $c/$o = 2L/^JT^l,l < L 

(I 7rL/4 if L <C l)-^^ We assume that the dot is 
not perturbed internally, so that the scattering is charac- 
terized by the energy-dependent scattering matrix S{e). 
Statistical averages over ensemble, denoted by (...), are 
found to leading order in 1 /N by diagrammatic technique 
in 5(e). ^'''^'^ We consider spinless electrons and normalize 
conductance/resistance by VsG^ jh^ where z/s = 2 is spin- 
degeneracy, and later compute the effect of SOI. ■^^^^^ 

The Hall resistance is measured as a linear response 
to additional voltages Vq applied at small frequency f2 
(experimentally, the frequency of lock-in < 100 Hz). 
If 1 and 2 are the source and drain, and the voltages 
are measured at 3 and 4, one has Ii — = I and 
/3 = /4 = deflning Rh = -Ri2,34 = (V3 - Vi)/h 
(Ref. 10) (applied voltages and current are related via 
Rccx = (Vi - V2)//i « (1/ni + l/n2)/N, whose fluc- 
tuations might be neglected here). A transposition (ex- 
change) of current- and voltage-probes, while keeping the 
same ac perturbations at to and the applied magnetic 
flux, gives a different Rjj = -R34,i2 = (Vi — V2)//3. Due 
to randomness of the voltage drop V34(i2) across the dot, 
the resistance is more relevant than the conductance con- 
sidered before. 

Derivation. For a non-interacting dot one can use scat- 
tering states (or Jauho-Meir-Wingreen formula '^°) to 
express the current in a contact a as a function of pertur- 
bations in a probe f3. It responds not only linearly to V/30, 
but also to Va^, and for small eKj/max {hu),T, ft/rd} <C 
1 we expand current to first orders in this small param- 
eter. This yields a conductance matrix gaf3 at ^ 0, 

8T r 

= g^^~- j '^^^^ [tc.S{e)tpS\e) - t^tp] fie) 
^ (eVp^V f, f{e + hu) + f{e-hcj)-2f{e) 

+ [-^) J (KJr 

xdetr laS{e)lpS\e) = ga0 + g'^p{eVi3^/2)\ (3) 

While the first term is the usual dc conductance g, the 
second is an out-of-equilibrium contribution. It is similar, 
but not identical, to the photo-assisted current. How- 
ever, Eq. (3) does not satisfy gauge-invariance: currents 
depend on V^^^ and change if all voltages are shifted by 
an arbitrary SV^^. Therefore, we consider the nearby gate 
and self-consistently find the internal potential of the dot 
U{t) due to linear screening: a potential Uf on some 
frequency / is a linear combination of voltages at this 
frequency, Uf = '^pupjVpj., where upj are complex 
characteristic potentials, which sum up to 1. This poten- 
tial depends on Td and RC time trc = ta/ {l + Vge^ /C/S). 
22,31 fpj^g overall voltage shift by — [/ reduces the problem 
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to the non-interacting, and now conductance, expressed 

-Ufi,Q)F^. (4) 



via = e^\V^ — C/p/4, reads as 



9at3 



9aP 



Since the sum of gaf3 in Eq. (4) over a or /3 vanishes, 
this sample-specific degenerate g satisfies charge conser- 
vation and gauge invariance, respectively: unlike Eq. (3), 
it depends on differences Vj — Vk unaffected by a voltage 
shift. Importantly, elements of g and their correlators 
depend both explicitly on the static Ufsfl and implicitly 
on dynamic u/j^i^. The matrix 52x2 oc 1 — is always 
symmetric, '^^ but the TRS breaking becomes noticeable 
in a multi-terminal setup. Indeed, due to the symmetry 
to the matrix transposition, g{(j)) — g^ {—(j)), the matrices 
f/(0) and (/"(O) are symmetric, but 34x4 is not. However, 
the measured result depends on the probe configuration. 

Indeed, if some i-th probe does not draw any current, 
we can eliminate the i-ih row and column from g and 
simultaneously shift all voltages by —Vi- We obtain Rh 
inverting the remaining part 33x3, and the same method 
gives Bifj in the other, transposed, configuration. When 
all V^i^ = 0, the substitute g ^ g reproduces 



p _ g3ig42 - g32g4i px _ gi3g24 - gi4g23 

nH,dc — -7— , ^H.dc — TT • 

det g3x3 det g3> 



3x3 



(5) 



and the symmetry of g immediately gives the Onsager re- 
lation, i?H,dc(0) = dc(~<^)- ^° These Hall resistances, 
always equal at (/) = 0, start to differ at weak magnetic 
flux 0^1, and to quantify their asymmetry we consider 
A = {Rh. dc — i?J^j.)/2(/). Its average vanishes, {A) — 0, 
and the Gaussian fluctuations are given by 



Var A = (A^) ~ (A)' 



{ni + n2)(rt3 -|- 714) 



IdT,(6) 



where X= (4/Td)e"^/^''(7rrr/n)2 sinh"^(7rrr/n). Com- 
pared to the classical Hall effect in such dots, this sen- 
sitivity is ~ {L/N\-pY ^ 1 stronger, but depends on 
the widths of the probes. As expected, when the dot is 
widely opened, <I>c $0 and Eq. (6) corresponds to the 
result of Ref. 28 up to a numerical coefficient ^ 1. 

With additional ac perturbations the Hall resistances 
[Eq. (5)] gain IFE-corrections, and we find that 5R± = 
SRh ± ^RJ{ are also normally distributed around zero, 
with 



Var SR± = 



(1-cosc.t)2 
Idr — — 



a/3 



(ni -l-n2)(n3 -1-714) 
x(Xi2+X34±2e-4^'^/"^ri2y34), (7) 



Interestingly, the aforementioned dependence of g and 
its fiuctuations on up Q vanishes from Eq. (7) due to the 
antisymmetry of Eq. (5). Therefore, the statistics of 5R± 
is unaffected by the static characteristic potentials. 
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FIG. 2: The width a (o) of the IFE Lorentzian distribution 
and (Tm (A) as functions of spin-orbit strength A, normalized 
by G of Eq. (8); (inset) rms of Rh (o) [A (□)], normalized 
by 8[16]/A''^, compared with theoretical dotted [solid] curves. 
Empty (filled) symbols correspond to A*' = 16(32). 



The magnetization IFE flux can now be evaluated us- 
ing the TRS breaking at = 0, <y9 = 5R^/2A. To 
leading order, we can take A and 5R- uncorrclatcd and 
using Eqs. (6) and (7) find a Lorentzian mesoscopic dis- 
tribution P{(p) — a I\k(ip^ -|- a'^)\ with 



/ dTl{l — coswr)^ 
4{huj)^ J drl 



[^12+^34-2ri2r34]3'e- («) 



The odd moments of P{ip) vanish and the even ones di- 
verge, but ip is solely defined by a, (\(p/a\^^^^) = \/2. 
In case one energy s among the energy scales fku, h/r^i, T 
is large compared to the other two, the ratio in Eq. (8) 
equals to l/e** multiplied by 3/8, 3/2 or 1/112 respec- 
tively, which results in Eq. (2). The cr^, a positive- 
semidefinite form of F-y, can be created by a single V 
in some contact and vanishes only in a degenerate sit- 
uation when all F^ are the same. One such example is 
when all voltages in the contacts have the same magni- 
tude and phase and a gate- voltage Vq ^ cannot induce 
IFE, see discussion of experiment. This situation is sim- 
ilar to a linearly polarized E not being able to induce 
the classical IFE, M cx [F x E*] = 0. Unless this uni- 
form regime is chosen to diminish Eq. (2) remains a 
good order-of-magnitude estimate of IFE. As a function 
of frequency, a in Eq. (8) is modified by screening: in a 
strongly interacting dot, trc <SC Td, in the high-frequency 
limit, ujtyic 3> 1, the capacitor is short-cut compared to 
the contact resistances and U k, Vq. As a result, IFE 
can occasionally become stronger than l/w^ due to the 
increased magnitude of — U\ for the particular con- 
figuration of voltages [see Fig. 1(b)]. 

Until now, spin-orbit was neglected in the derivation 
of Eq. (8) . To take it into account we construct a scat- 
tering matrix S{e), which depends on magnetic field and 
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SOI-strcngth, compute g and g, and find IFE flux (p = 
6R-/2A. For illustration, wc take a simplified model 
of SOI using a parameter A ^ H/Atso (Ref. 26) for the 
dot's random Hamiltonian ?i: in the limits A = 0(oo) it 
belongs to the Gaussian orthogonal (symplectic) Ensem- 
ble GOE, P = I (GSE, /? = 4). In the crossover region 
the 2M X 2M Hamiltonian with a fixed A is represented 
hy n = [VMHi + (A/2)H4]/VmTA2. Our numerics 
are done at M = 25N to ensure M > TV, for T = 0, 
w < h/Td, Ni = N/4, and only Vi ^ 0. To find the 
IFE scale we fit 1200 sample statistics to a shifted 
Lorentzian, and the result is presented in Fig. 2 together 
with a given by moments, am — ily^l^^^) / {\'p\~^^'^) ■ De- 
viations of (7(A) from Eq. (8) are attributed to a rela- 
tively small number of channels and appear to be non- 
systematic. The inset in Fig. 2 shows rms of Rh and A 
from the fits to normal distributions in very good agree- 
ment with our predictions 0.5(l-h3/(l-h4A7iV)2)i/2 and 
0.5(l + (3-4AViV)/(l-h4AV7V)4)i/2, respectively. The 
lack of any systematic trend of the available statistics in 
the main plot (compared to the inset) suggests that the 
mesoscopic IFE is unaffected by SOI, in contrast with 
the weak IFE existing only due to SOI. ^ 

The generalization from a 4-contact setup to M con- 
tacts includes arbitrarily coupled floating contacts, which 
are typical in a Hall measurement. The voltage in a 
floating probe is adjusted to allow exchange of electrons 
with the dot, but not to draw any current. As a re- 
sult, our analytical results become renormalized using 
A'' = J2t=i + ^f- -^ot necessarily integer, Nf is de- 
fined similarly to Ni by the total dimcnsionlcss conduc- 
tance Nf of all available fioating contacts. In the dephas- 
ing probe model the inelastic scattering with time 
adds h/ATcf, fictitious channels to Nf. This redefinition 
decreases Td and Ua used in 2 and Eqs. (6) (8). 

Discussion of experiment. Reference 14 measured a 
TRS breaking at </> = in gated ballistic Hall sam- 
ples, where each of 6 (=4-|-2 floating) contacts had 
Na ~ 200 channels. AC perturbations were applied at 
CO lO*' — 10^°Hz, either asymmetrically (sample A) or 
uniformly (sample B). Corrections to Rh, quadratic in 
small perturbation amplitude, were measured as func- 
tions of w in the sample A. They vanish for w < wq and for 
CO > u)c become of the same order, and a single-parameter 
pumping in a billiard was used to numerically reproduce 



non-monotonic Rh — R"h as a function of ui. On the con- 
trast, sample B showed no TRS breaking, 6Rh ~ SRjj. 

Alternatively, data can be interpreted as mesoscopic 
fluctuations of IFE in a chaotic sample. They are visible 
even at large N due to vanishing classical effect. Unfor- 
tunately, data can not be directly compared with RMT 
result (7) because the Hall cross docs allow direct trajec- 
tories resulting in much smaller non-local voltages and at 
0.3 K it is beyond the universal regime. Experiment inter- 
prets LOo as a frequency characteristic for (unknown) ca- 
pacitive coupling with the contacts, when ac voltages be- 
come noticeable. In our setup the voltages are given and 
this circuit effect is not accounted for. At low frequencies 
5R± vary cx (uJTdY'^, but the anti-symmetric component 
is generally smaller by construction [cf. ± in Eq. (7)]. 
Beyond the threshold lOc ^ I/t^ both fluctuate similarly 
with a typical period Ato ~ l/ra. Indeed, in the exper- 
iment is usually larger and lUc ^ Alu is consistent 
with Td «3 ps expected from a dot with the sample area 
and N w 1200. For the uniformly perturbed sample B 
with equal voltages V1...4, Eq. (7) indeed results in fluc- 
tuations of (5_R-|_ and SR- = observed in experiment, 
which does not enter a high frequency regime. How- 
ever, it is desirable to perform measurement in a truly 
chaotic sample without direct trajectories and at lower 
T and N to increase quantum fluctuations, or measure 
IFE directly using superconducting quantum interference 
device, similarly to persistent current. 

Conclusions. We consider a time-reversal symmetry 
breaking by external ac voltages. Rectified perturbations 
generate a static magnetoelectric effect, a random mag- 
netization of a mesoscopic sample, similar to the classi- 
cal Inverse Faraday Effect. Mesoscopic fluctuations of 
the magnetization flux can be measured using out-of- 
equilibrium transport in a multi-terminal quantum dot. 
We estimate a typical flux, quadratic in voltag 
function of frequency, screening, coupling to reservoirs, 
and numerically find that spin-orbit scattering has a very 
small effect on IFE. Our results allow qualitative compar- 
ison with and explain most important features of a recent 
experiment in TRS breaking. 
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